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Abstract 

We search for the decays 5° ^ p°p°^ B° /9°/o(980), and B° /o(980)/o(980) in a sam- 
ple of about 348 milhon r{4S) BB decays cohected with the BABAR detector at the PEP-II 
asymmetric-energy e~^e~ collider at SLAC. We find evidence for B^ p^p^ with 3.0a significance 
and measure the branching fraction B = {1.16~^q'^q zb 0.27) x 10~^ and longitudinal polarization 
fraction fi = OMtlH ± 0.05. As a consequence, the uncertainty on the CKM unitarity angle a 
due to penguin contributions in B ^ pp decays is estimated to be 18° at la level. We also set upper 
limits on the B^ — > p^fo{980) and B^ — > /o(980)/o(980) decay rates. All results are preliminary. 
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1 INTRODUCTION 



Measurements of CP-violating asymmetries in the B^B^ system provide tests of the Standard Model 
by over-constraining the Cabibbo-Kobayashi-Maskawa (CKM) quark-mixing matrix HI through the 
measurement of the unitarity angles. Measuring the time-dependent CP asymmetry in a neutral- 
5-meson decay to a CP eigenstate dominated by the tree-level amplitude b uud gives an approx- 
imation aeflf to the CKM unitarity angle a = arg [—V^^V^l/V^^V*jJ\. The correction Aa = a — Oes 
accounts for the additional contributions from loop (penguin) amplitudes. The value of Aa can be 
extracted from an analysis of the branching fractions of the B decays into the full set of isospin- 
related channels [^j. 

Measurements of branching fractions and time-dependent CP asymmetries in B ^ vrvr, pn, 
and pp have already provided information on a. Because the branching fraction^ for B^ tt^tt^ 
is comparable to that for B~^ — > 7r~^Tr^ and B^ — > 7r"^7r~, the limit on the correction is weak: 
|Aa7r7r| < 41° at 90% confidence level (C.L.) On the contrary, the B^ —>■ p^p^ decay has a much 
smaller branching fraction than B^ — > p'^p~ and B~^ P^P^ channels |ll-llOj. As a consequence, it 
is possible to set a tighter limit on Aopp [21IZ1^B- This makes the pp system particularly effective 
for measuring a. 

In B ^ pp decays the final state is a superposition of CP-odd and CP-even states. An isospin- 
triangle relation [2] holds for each of the three helicity amplitudes, which can be separated through 
an angular analysis. The measured polarizations in B~^ P^P^ IS lEl and B^ — > p'^p~ [ZHHl 
modes indicate that the p's are nearly entirely longitudinally polarized. In this paper we present 
evidence for the B^ p^p^ decay, the first measurement of the longitudinal polarization fraction 
in this decay, and updated constraints on the penguin contribution to the measurement of the 
unitarity angle a. These results supersede our previous limits on this decay [HISl- 

2 THE BABAR DETECTOR AND DATASET 

These results are based on data collected with the BABAR detector |12| at the PEP-II asymmetric- 
energy e~^e~ collider ^H] located at the Stanford Linear Accelerator Center. A sample of 347.5 ± 
1.9 million BB pairs, corresponding to an integrated luminosity of approximately 316 fb~^, was 
recorded at the T(45) resonance with the center-of-mass (cm.) energy ^/s = 10.58 GeV. We use 
a sample of 28 fb~^ taken 40 MeV below the T{AS) resonance to study background contributions 
from e~^e~ qq (q = u, d, s,or c) continuum events. Charged-particle momenta and trajectories 
are measured in a tracking system consisting of a five-layer double-sided silicon vertex tracker 
and a 40-layer drift chamber, both within a 1.5-T solenoidal magnetic field. Charged-particle 
identification is provided by measurements of the energy loss in the tracking devices and by a 
ring- imaging Cherenkov detector. 

3 ANALYSIS METHOD 

We select B^ — > p^p^ (7r+7r~)(7r+7r~) candidates from combinations of four charged tracks 
that are consistent with originating from a single vertex near the e~^e~ interaction point. The 
identification of signal B candidates is based on several kinematic variables. The beam-energy- 
substituted mass, mES = [(•s/2 + Pj •ps)^/-Ef — p^]"^''^, where the initial total e+e~ four-momentum 

^Charge conjugate B decay modes are implied in this paper. 
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{Ei, Pi) and the B momentum pB are defined in the laboratory frame, is centered near B mass 
with a resolution of 2.6 MeV for signal candidates. The difference between the reconstructed B 
energy in the cm. frame and its known value AE = E"^^ — \/s/2 has a maximum near zero with a 
resolution of 20 MeV for signal events. Four other kinematic variables describe two possible tt~^tt~ 
pairs: they are invariant masses mi(7r7r) = mi and m2{'KTi) = m2, and helicity angles 9i, 02- 

The angular distribution of the B^ p^p^ decay products can be expressed as a function 
of the helicity angles 9i and 02, which are defined as the angles between the direction of vr"'" 
and the direction of the B in the rest system of each p^. The resulting angular distribution 
(PT/{T dcosOi dcos6i2) is 

^ - /l) sin^ 01 sin^ 02 + h cos^ Oi cos^ 02 1 , (1) 

where fi = \Aq\'^/{'E\Ax\'^) is the longitudinal polarization fraction and Ax=-ifl^^i are the helicity 
amplitudes. 

The selection requirements for signal candidates are the following: 5.240 < rriEs < 5.290 GeV/c^, 
lAi^l < 85 MeV, 550 < mi^2 < 1050 MeV, and | cos0i^2| < 0.98. The latter requirement removes a 
region with low reconstruction efficiency. In addition, we veto the copious decays B^ D^*^^7r^ — > 
{h~^TT~7r~)7r~^ , where refers to a pion or kaon, by requiring the invariant mass of the three-particle 
combination to differ from the D-meson mass by more than 13.2 MeV, or 40 MeV if the kaon is 
positively identified. 

We reject the dominant continuum background by requiring | cos 0^1 < 0.8, where 6t is the 
angle between the S-candidate thrust axis and that of the remaining tracks and neutral clusters 
in the event, calculated in the cm. frame. We suppress continuum background further using the 
polar angles of the B momentum vector and the S-candidate thrust axis with respect to the beam 
axis in the cm. frame. Other discriminating variables calculated in the cm. frame include the two 
Legendre moments Lq and L2 of the energy flow around the -B-candidate thrust axis 14^ and the 
sum of the transverse momenta of all particles in the rest of the event, calculated with respect to 
the B direction. These variables are combined in a neural network, whose output is transformed 
into the approximately Gaussian-distributed variable S. 

After application of all selection criteria, A^cand = 65180 events are retained, most of which are 
background events, well separated in the kinematic observables from B^ — > p^p^, B^ — > p^fQ{980), 
and B° /o (980) /o (980) signal candidates. On average, each selected background event has 1.05 
candidates, while in Monte Carlo (MC) samples we find 1.15 and 1.03 candidates for longitudinally 
and transversely polarized B^ — > p^p^ decays, respectively. When more than one candidate is 
present in the same event, the candidate having the best consistency with a single four-pion 
vertex is selected. 

The signal selection efficiency determined from Monte Carlo J^j simulation is 23.5% or 28.9% 
for longitudinally or transversely polarized events, respectively. MC simulation shows that 18% 
of longitudinally and 4% of transversely polarized signal events are misreconstructed with one or 
more tracks not originating from the B^ — > p^p^ decay. These are mostly due to combinatorial 
background from low-momentum tracks from the other B meson in the event. 

Further background separation is achieved by the use of multivariate i?-flavor-tagging algo- 
rithms trained to identify primary leptons, kaons, soft pions, and high-momentum charged particles 
from the other B ^j. The discrimination power arises from the difference between the tagging 
efficiencies for signal and background in seven tagging categories (ctag = 1 ~ 7). 
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4 MAXIMUM LIKELIHOOD FIT 



We use an unbinned extended maximum likelihood fit to extract the p^p^ event yield and 

fraction of longitudinal polarization f^. We also fit for the event yield of p^fo and /q/q 

decays, as well as several background categories. The likelihood function is 



where Uj is the number of events for each hypothesis j (signal B^ p^p^ , five other i3-decay classes, 
and continuum), and Vj{xi) is the corresponding probability density function (PDF), evaluated with 
the variables Xi = {mES; ^-^'j "ii, m2, cos ^i, cos 6*2, Ctag} of the ith event. 

We use MC-simulated events to parameterize contributions from other B decays. The charmless 
modes are grouped into several classes with similar kinematic and topological properties: B^ — 
^"70(990); B'^ — > /o(980)/o(980); B'^ — > afir^; and a combination of other charmless modes, 
including B^ — > p^K*'^, B^ P^P^^ B pn, and B^ — > p^p^- One additional class accounts for 
the remaining neutral and charged B decays to charm modes. The number of events in each class 
Uj is left free in the fit. We ignore any other four-pion final states in our invariant mass window 
whose contributions are expected to be small. 

Since the statistical correlations among the variables are found to be small, we take each Vj 
as the product of the PDFs for the separate variables. Exceptions are the kinematic correlation 
between the two helicity angles in signal, and mass-helicity correlations in other 5-decay classes 
and misreconstructed signal. They are taken into account as discussed below. 

We use double-Gaussian functions to parameterize the mss and A.E PDFs for signal, and a 
relativistic Breit-Wigner (BW) for the resonance masses of p^ jlH] and /o(980) ^Hl- The angular 
distribution at production for B^ — > p^p^, B^ P^fo, and B^ — > /o/o modes (expressed as a 
function of the longitudinal polarization in Eq. for B^ P^P^) is multiplied by a detector 
acceptance function ^(cos ^i, cos 6*2), determined from MC. The distributions of misreconstructed 
signal events are parameterized with empirical shapes in a way similar to that used for B back- 
ground, as described below. The £ variable is described by three asymmetric Gaussian functions 
with different parameters for signal and background distributions. 

The PDFs for exclusive non-signal B decay modes are generally modeled with empirical ana- 
lytical distributions. Several variables have distributions identical to those for signal, such as rriEs 
when all four tracks come from the same B, or 7r~''7r~ invariant mass mi^2 when both tracks come 
from a p^ meson. In certain exclusive modes the two candidates can have very different mass 
and helicity distributions, e.g. when only one of the two p^ candidates is a genuine meson or 
when one of the two p^ candidates contains a high-momentum pion (as in — > aivr). In such cases, 
we use a four- variable correlated mass-helicity PDF. 

The signal and i?-background PDF parameters are extracted from MC simulation. The initial 
continuum background PDF parameters are obtained from data in mss and AE sidebands and are 
then left free in the fit. The MC parameters for ttiesi AE, and £ PDFs are adjusted by comparing 
data and MC in control channels with similar kinematics and topology, such as B^ D~7r~^ with 
D~ — > K'^TT^TT^ . Finally, the i?-flavor tagging PDFs for all decay modes are the normalized discrete 
Ctag distributions of tagging categories. Large samples of fully reconstructed i?-meson decays are 
used to obtain the -B-tagging efficiencies for signal B decays and to study systematic uncertainties 
in the MC values of i3-tagging efficiencies for the B backgrounds. 
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5 RESULTS 



Table ^ shows the results of the fit. The — > p^p^ decay is observed with a significance of 
3. Oct, as determined by the quantity \/—2 log(£o/>Cmax), where £max is the maximum likelihood 
value, and Co is the likelihood for a fit with the signal contribution set to zero. It corresponds to 
a probability of background fluctuation to the observed signal yield of 0.1%, including systematic 
uncertainties, which are assumed to be Gaussian-distributed. We do not observe significant event 
yields for 5° ^ p°/o(980) and 5° ^ /o(980)/o(980) decays. Background yields are found to be 
consistent with expectations. In Fig. ^ we show the projections of the fit results onto tit-es and AE. 

Table 1: Summary of results: signal yield (ngig, events), fraction of longitudinal polarization (/l), 
selection efficiency (Eff), branching fraction (^sig), branching fraction upper limit (UL) at 90% CL, 
and significance (including systematic uncertainties). The systematic errors are quoted last. We 
also show the background yields for aivr, BB, and qq components (events, with only statistical 
uncertainties quoted). 



Quantity 




Value 


nsig (5° - pV°) 
/l 




98t^i ± 22 
0.86l|j:JJ ± 0.05 


Eff (%) 




24.2 ± 1.0 


Bsig (xlO 




l-16Io.36 ±0-27 


Significance (a) 




3.0 (3.4 statistics only) 


nsig (5° ^ pO/o) 




Utlj ± 13 


Eff (%) 




20.7 ±0.8 


^sig X B{fo vr+vr" 


) (xlO-6) 


0.17t[J:i ±0.18 


ULxS(/0 ^ TT+TT-) 


(xlO-6) 


0.68 


r^sig (i3° ^ Mo) 




-5tl ± 12 


Eff (%) 




23.5 ±0.9 


Bsig X B^ifo vr+vr 


-) (xlO-6) 


-0.06l[j:|}^ ±0.15 


ULxfi2(/o ^ vr+^-; 


) (xlO-6) 


0.33 


'^charmless 
^99 




9011 

17+113 
"-■■'-99 

328011^7 
6171911^ 
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niEs'GeV) AE(GeV) 



Figure 1: Projections of the multidimensional fit onto m-Es and AE after a requirement on the 
signal-to-background probability ratio with the plotted variable excluded. This requirement maxi- 
mizes the fraction of signal events in the sample. The data points are overlaid by the solid (dashed) 
line, which shows the full (background only) PDF projection. The individual PDF components are 
shown for the —>■ p^p^ (solid red) and —>■ p^fo modes (dotted green). 

6 SYSTEMATIC STUDIES 

Dominant systematic errors in the fit originate from statistical uncertainties in the PDF parame- 
terizations, due to the limited number of events in the control samples. The PDF parameters are 
varied by their respective uncertainties to derive the corresponding systematic errors (15, 11, 12 
events for p^p^, p^fo, and /o/o respectively, and 0.05 for f^). We also assign a systematic error 
of 2 events for p^p^ and /o/o and 7 events for p^fo (0.02 for f^) to account for a possible fit bias, 
evaluated with MC experiments. The above systematic uncertainties do not scale with event yield 
and are included in the calculation of the significance of the result. We also assign 8%, 5%, 10% 
multiplicative systematic error due to possible fit bias for p^p^, p^fo, and /o/o modes, respectively. 

We estimate the systematic uncertainty due to the interference between p'^p^ and afn^ final 
states using simulated samples in which the decay amplitudes for B^ — > p^p^ are generated ac- 
cording to this measurement and those for B^ — > a^vr^ correspond to a branching fraction of 
(39.7 ±3.7) X lO-'^ [10]. Their amplitudes are modeled with a Breit-Wi gner function for all p — > tttt 
and ai pn combinations and their relative phase is assumed to be constant across the phase 
space. The strong phases and CP content of the interfering state afn^ are varied between zero and 
a maximum using uniform prior distributions. We take the RMS variation of the average signal 
yield (14 events for the p^p^ yield, or 0.03 for f^) as a systematic uncertainty. 

Uncertainties in the reconstruction efficiency arise from track finding (2%), particle identification 
(2%), and other selection requirements, such as vertex probability (2%), track multiplicity (1%), 
and thrust angle (1%). 

7 IMPLICATIONS FOR THE CKM ANGLE a 

Since the tree contribution to the B^ —>■ p^p^ decay is color-suppressed, the decay rate is sensitive 
to the penguin amplitude. Thus, this mode has important implications for constraining the uncer- 
tainty in the measurement of the CKM unitarity angle a due to penguin contributions to i? — > 
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Figure 2: Ax^ on Aopp obtained from the isospin analysis discussed in the text. The dashed Unes 
at Ax^ = 1 and Ax^ = 2.7 are taken for the la (68%) and 1.64(j (90%) interval estimates. 



decays. 

In the isospin analysis 0) minimize a x^ that includes the measured quantities expressed 
as the lengths of the sides of the isospin triangles. We use the measured branching fractions and 
fractions of longitudinal polarization of the P^P^ [El ^1 and P^P~ [HI HI decays, 

the CP-violating parameters and obtained from the time evolution of the longitudinally 
polarized B^ — > p'^p~ decay 011], and the branching fraction of B^ — > p^p^ from this analysis. 
We assume Gaussian behavior of the distributions. We neglect 1 = 1 isospin contributions, non- 
resonant and isospin-breaking effects. 

With the — > p'^p^ measurement we obtain the constraint on a due to the penguin contribu- 
tion and obtain a 68% (90%) CL hmit on Aopp = a - a^a of ±18° (±21°). Fig. Ulshows the Ax^ 
on Aopp. The central value of a obtained from the isospin analysis is the same as Oes, which is 
constrained by the relation sin(2aefr) = /(^ ~ C'^~^)^/^ and is measured with the B'^ P^ P~ 
decay 

The error due to the penguin contribution becomes the dominant uncertainty in the measure- 
ment of a using B ^ pp decays. However, once the sample of B^ — > p^p^ decays becomes more 
significant, time-dependent angular analysis will allow us to measure the CP parameters and 
Cj^, analogous to 5"^" and , resolving ambiguities inherent to isospin triangle orientations. 
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8 SUMMARY 



In summary, we have found evidence for p^p" decay with 3.0a significance. We measure 

B{B^ p^p°) = l.ietUI^ (stat.) ±0.27 (syst.) 

and we determine the longitudinal polarization fraction for these decays of 

= OMtoil (stat.) ± 0.05 (syst.) 

The measurement of this branching fraction combined with that for P^P^ smd B^ P^P~ 

decays provides a constraint on the penguin uncertainty in the determination of the CKM unitarity 
angle a. We find no significant evidence for the decays B^ — > p^fo and B^ — > /o/o- These results 
are preliminary, and they supersede our previous measurements [3] . 
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